Evolution of pairwise entanglement in a coupled n— body system 
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We study the exact evolution of two non-interacting qubits, initially in a Bell state, in the presence 
of an environment, modeled by a kicked Ising spin chain. Dynamics of this model range from 
integrable to chaotic and we can handle numèrics for a large number of qubits. We find that the 
entanglement (as measured by concurrence) of the two qubits has a close relation to the purity of 
the pair, and follows an analytic relation derived for Werner states. As a col·lateral result we find 
that an integrable environment causes quadratic decay of concurrence, while a chaotic environment 
causes linear decay. 



Ever since the famous paper of Einstein, Podolski, and 
Rosen entanglement and Bell states are at the cen- 
ter of interest of those concerned with the foundations 
of quantum mechanics. Recent experimental realizations 
of deterministic teleportation with photons 0, and 
ions 0, 0] and a rising technical interest in quantum 
information 0| have stirred a great deal of interest in 
the decay of pairwise entanglement due to decoherence. 
Decoherence is one of the main problems in such appli- 
cations 0, H, 01 j tnc other being systematic errors in the 
physical implementation of a logical or mathematical al- 
gorithm [ïfl[ïlE3· 

In the present letter we shall present a model involv- 
ing many qubits, that is well suited to study this kind of 
qüestions, and then proceed to apply it to the evolution 
of an initial Bell state multiplied by a random state in 
cnvironments with different dynamical properties. This 
random state can be interpreted as a high temperature 
environment or as reflecting our ignorance as to the state 
of the environment. We shall consider both the behav- 
ior of purity and concurrence of the pair under this evo- 
lution, and determine the influence of dynamics of the 
environment, in particular the question whether integra- 
bility or chaoticity of the environment are relevant. We 
shall find, that this is indeed the case, but beyond that 
we shall find, that the average relation between purity 
and concurrence is independent of the dynamics. Thcrc 
is also strong indication, that self averaging will tend to 
reduce the variance of the ensemble, as we increase the 
number of qubits that rcprcscnt the environment. 

The model itself is a generalization of the kicked spin 
chain introduced by Prosen [Rg and used in a similar 
context as the present one to study decoherence or en- 
tanglement under echo-dynamics The model had 
the virtue of allowing numèrics with rather large num- 
bers of qubits and in integrable situations some analytic 
solutions. The original model left no frecdom to choose 
the coupling strength between different parts of the sys- 
tem. The ease of numerical handling of the chain does 
not depend upon this fact. We can thus generalize to 
arbitrary two-body couplings. This allows us to model 



different situations. We can thus consider a central sys- 
tem consisting of two spins weakly coupled to one or two 
spin chains acting as environment, whose dynamics may 
vary from integrable to chaotic [,15]. Because we want 
the concurrence to be affected exclusively by the cou- 
pling to the environment we choose non-interacting spins 
for the selected pair. For this environment it is sensible 
to consider random states to emulate a bath at fairly 
high temperature. Using unitary time evolution of the 
total system and partial tracing over the environment we 
can then calculate concurrence and purity decay of the 
selected pair, and discuss their behavior. Their depen- 
dence on the dynamics of the environment is significant. 
Yet the relation between purity and concurrence decay 
that is known for Werner states will be seen to hold to 
good approximation in all dynamics studied. 

The concurrence (C) can be regarded as a good mea- 
sure of entanglement for a density matrix p 16]. Con- 
currence is defined as 

C = max{0, Ai - A 2 - A 3 - A 4 } (1) 

where Ai are the eigenvalues of the matrix 
\J p{p y ® Vy)P*(?y ® o'y) m non-increasing order, 
(*) denotes complex conjugation in the computational 
basis and a y is a Pauli matrix. 

Purity is a measure of mixedness for a state character- 
ized by a density matrix p in an arbitrary Hilbert space. 
It has a value of one for pure states and less than one in 
any other casc. It is defined as 

P(P) = Trp 2 (2) 

We use this measure instead of the usual von Neumann 
entropy (defined as — Trplogp) since it is easier to handle 
from an algebraic point of view and both measures, albeit 
different, contain similar information. 

Having a given value of purity in general does not fix 
the value concurrence and vice versa. We can visualizc 
the set of physical states in a plane with C and P in the 
axes. In Fig. Elthe gray region plus the line {(P,0)\P € 
[1/4, 1/3]} represent the range of compatible vàlues for 
concurrence and purity |l7j . 



Let us define a general Bell state as I^Bcn) = 
(Iaíi) |M2> + \m)\m)) I V% where {\Hi),\r]i)} define an or- 
thonormal basis for particle i. We shall now proceed 
to calculate the relation between purity and concurrcncc 
for a Werner state in which the entangled component is 
a general Bell state: 
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where a lies between zero and one. 

Taking into account that purity and concurrence do 
not change under any independent particle unitary trans- 
formation, we may choose the computational basis so as 
to write |V>Bdi) = (100) + |H))/\/2 in Eq. yielding the 
cxplicit form of the density matrix 
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Then we calculate the concurrence (Eq. ^| and purity 
(Eq. [21, obtaining the exact expressions P = 1 — 3a/2 + 
3a 2 /4 and C — max{0, 1 — 3a/2}. Thus concurrence is 
given in terms of the purity by 
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if 1/3 < P < 1 
if 1/4 < P < 1/3. 



(5) 



The Hamiltonian of the kicked Ising model is 
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h = jj,k<Tj(Tk 

j,k=0 

(6) 

where <5i (t) represents an infinite train of Dirac delta 
functions with period one. The model thus consists of 
a periòdic chain of spin 1/2 partides coupled to all other 
spins by an Ising intcraction (first term) and periodi- 
cally kicked by a tilted magnètic field (second term). 
Due to the 2 body nature of the Hamiltonian we are 
able to calculate the time evolution of arbitrary ini- 
tial conditions for up to 24 qubits. Setting = 1.4, 
Jj.k = Sj+i,k and varying the paral·lel component of the 
magnètic field we can obtain integrable (6|| = 0), non- 
ergodic and non-integrable (0 < b\\ S 0.8), and fully er- 
godic (0.8 ^ ò|| íí 1.4) dynamics |13| . Our model differs 
from the one given in |l3| by the fact that the coupling 
Jj k is between any pair of partides and has arbitrary 
strength, instead of Jj t k — JSj+i,k, which couples near- 
est neighbors with fix strength. Throughout this letter 
we will couple only first neighbors, but we shall keep the 
freedom of choosing arbitrary strength. 

To considcr the situation described in the beginning of 
the letter we must weakly couple 2 spins, say spins "0" 
and "1" , to the rest of the chain, which we will consider as 
the environment i. e. J\ and Jl-i much smaller than the 
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FIG. 1: Different configurations of the coupling of the system 
to the environment. The open circles represent the central 
system and the filled circles th environment. Thick/thin lines 
represent strong/weak interaction. In (a) we see both par- 
tides coupled to one environment, in (b) only one particle 
coupled to the environment and in (c) each particle coupled 
to two independent environments. 



typical Ising interaction within the environment, which 
we choose to be 1. We also set Jo = in order to prevent 
any interaction between the spins in the central system. 
The fact that we keep the kick in the central system can 
represent local operations made by the "owners" of each 
of the qubits, and will not affect the vàlues of concurrence 
and purity. 

Rewriting the Hamiltonian in terms of central system, 
environment and interaction as Pi = H c + H c + H ce the 
parts are given by 



H c = St(t) J>xo? + &||0?), 



(7) 
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H e = J2 Jrf°3+i + Si (t) ( 6 W + ò ll (8) 

J=2 j=2 

H cc = Jl-i(t z L -iVq + Jicrlal- (9) 

We shall consider three particular situations. The first 
one represents both partides coupled with equal strength 
to one environment, i. e. Jo = 0, Jl-i — Ji an d 
Ji = ■■■ = Jl-2 = 1) see Fig. nja). In the second 
one only one particle is coupled to the environment i. e. 
Ja = Jl-i = and J 2 = • • • = Jl-2 = 1, see Fig. IHb). 
Finally we can couple each particle to independent envi- 
ronments setting, e. g., Jl-i = Ji, Jo = Jk = for some 
2 < k < L — 2 and all the other J's equal to one, see 
Fig. ^c). The non- unitary evolution of the central sys- 
tem is calculated performing a unitary evolution of the 
whole system (yielding state \ip(t)}) and then performing 
the partial trace over the environment, that is 



p(t) = Tr cnv \iP(t))m)\- 
We shall consider an initial condition 



|^(í = 0)) = I^Bell) ® \t 
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i. e. a general Bell state not entangled with the envi- 
ronment. Note that the reduced density matrix of the 
central system at this time is IV'BeiiXV'Beiil (a pure Bell 
state) which implies C(p(0)) = P(p(0)) = 1. 
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FIG. 2: (Color online) This figures shows the evolution 
of concurrence for 2 different thermal baths, with 60 initial 
condition given by Eq. 1111 The thermal bath is either chaotic 
(with a 6|| = b± — 1.4, green continuous lines) or regular 
(6ii = and b± — 1.55, red dashed lines). The corresponding 
points indicate the average over the initial conditions. We 
have 15 qubits and set JcoupUng = 0.03. 
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FIG. 3: (Color online) This figures shows the evolution 
of the system in the (C, P) plane, for different dynami- 
cal regimes: integrable (h± = 1.55, hu = 0), intermediate 
(h± « 1.89, h\\ » 0.59), and chaotic (h± = 1.4, h\\ = 1.4), in 
blue triangles, green squares, and orange circles respectively. 
In this figure L = 14 and J c = 0.01. We use 10 initial condi- 
tions and a total time of 4500 steps. The red line shows the 
relation for Werner states Eq.|S] 



We now present the results of our mimcrical calcula- 
tions, of both concurrence and purity of the selected pair 
of spins as a function of time. We first concentrate on 
configuration (a). 

Fig. |21 shows the time evolution of concurrence of the 
selected pair of spins for the initial state consisting of a 
Bell pair coupled to random stated for the environmcnt. 
We choose the environment in one case to be ergodic and 
in the other to be integrable using corresponding param- 
eters from Ref. ,13] as mentioned above. It is pertinent to 
mention, that these dynamical properties were obtained 
for a cyclic chain, while the chains representing our en- 
vironmcnts are open. Yet for large numbers of spins this 
is irrelevant. With the first configuration, i. e. an en- 
vironment of 13 spins coupled at each end to one of the 
selected spins, that form the central system; 60 states are 
chosen at random as in Eq. 1111 Each of the sixty mem- 
bers of the ensemble as well as the average is plotted. 
Note that for integrable environment all curves coincide 
within line width; it is not clear whether this difference 
in variance between the chaotic and the integrable case is 
genèric, and we plan further investigation on that point. 

Concurrence decays considerably faster for an inte- 
grable environment, than for an ergodic one. This fact 
is remarkable, as it is quite different from findings in 
Ref. . We observe a similar behavior for the purity of 
the central system; this is entirely in keeping with find- 
ings for echo dynamics with a similar kicked spin chain, 
though in this case echo dynamics with strong coupling 
between the central system and the environment was con- 
sidered. Indeed, just as for purity in echo dynamics, be- 
yond the Zenon time (around 2 time steps) the decay is 
linear in the ergodic case and quadratic in the integrable 
one. This is not surprising, as the linear response result 



derived in |14| for the decay of purity in echo dynamics 
trivially translates to purity in forward time evolution, if 
we consider the coupling, Eq. [ü] to be the perturbation 
of the echo. This is also in keeping with linear response 
calculations and numèrics for entanglement production 
given in Ref. |l9l . |2(| . The result is not in contradiction 
to contrary findings for coherent states [2ll for inte- 
grable systems as these have decay times governed by a 
different h dependence ■ 

While the results presented are by thcmselves of con- 
siderable interest as they are somewhat counterintuitivc, 
we now wish to focus on a different aspect, namely that 
a one to one relation between concurrence and purity 
emerges. We thus plot in Fig. [3] concurrence versus pu- 
rity for 14 qubits, averaging purity for a given concur- 
rence, again over 10 initial conditions of the form given 
in Eq. ^2 we pl°t the results both for the integrable, the 
ergodic and an intermediate choice of the environment. 
We choose only 14 spins because with a bigger system 
the curves would become indistinguishable in the figure. 
Remarkably we find that the three cases coincide within 
statistical error and that they agree with the relation 
given in Eq. [3] for Werner states, though they dennitely 
are not Werner states; the latter fact was checked directly 
by considering the eigenvalues of the density matrix. We 
thus find, that for a random environment on average the 
relation[S]holds though it was derived without dynamics. 

Due to self-averaging this seems also to oceur for an 
individual typical random state of the environment as the 
number of qubits increases; furthermore the range over 
which the relation holds also increases with the size of the 
environment, as we can see in Fig. ^Iwhere the relation 
is plotted for different numbers of qubits. 
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FIG. 4: (Color online) This figures shows the evolution for 
one initial condition in time steps of 125, in the (C, P) plane, 
in the chaotic regime (òii = b± — 1.4). The number of qubits 
is in the reservoir is varied, the red triangles, purple circles, 
blue squares, and green rhombus correspond to 5, 8, 11, and 
14 qubits. The picture suggests that increasing the size of 
the environment improves the Werner state approximation. 
Similar results are obtained in the regular regime. 

All results presented are calculated for the first of the 
three configurations shown in Fig. Q The results for the 
other configurations are quite similar, after proper rescal- 
ing of the coupling strength. We therefore do not show 
them. Nevertheless they are of interest. For the second 
configuration one partner of the initial Bell pair is dy- 
namically speaking a spectator, yet the evolution of con- 
currence and purity of the pair are marginally affected. 
The last case is more of an instructive toy model: Here 
we have two uncoupled environments, and we can start 
with a random state in each of these. The purity of the 
uncoupled subsystems will remain unchanged, but the 
purity and concurrence of the initial Bell pair will decay. 
Thus one might consider seeing a paradox, but this is not 
the case; the entanglement of the pair simply is spread 
over all the system with time. Though the three config- 
urations are physically quite different and the individual 
behavior of purity and concurrence is slightly affected, 
the relation between the two is entircly robust. 

Summarizing, we have coupled a non-interacting Bell 
pair to an environment, that allows dynamics to be varied 
from integrable to chaotic in a smooth way, in the frame- 
work of a generalized kicked Ising chain. This model 
allows fairly large calculations, and we have found that 
under a wide variety of couplings and dynamics for the 
environment the relation between the purity and the con- 
currence of the Bell pair along its time evolution follows 
quite closely the one known for Werner states. ff we 
average over random environment states and choose a 
sufficicntly large number of qubits we follow this relation 
almost exactly. As a side-product we found that over the 
entire range of environment dynamics concurrence for an 
integrable environment decays faster than for a chaotic 



one. 
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